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In the context of a microscopic model of string-inspired foam, in which foamy structures are 
provided by brany point-like defects (D-particles) in space-time, we discuss flavour mixing as a 
result of flavour non-preserving interactions of (low-energy) fermionic stringy matter excitations 
with the defects. Such interactions involve splitting and capture of the matter string state by the 
defect, and subsequent re-emission. As a result of charge conservation, only electrically neutral 
matter can interact with the D-particles. Quantum fluctuations of the D-particles induce a non- 
trivial space-time background; in some circumstances this could be akin to a cosmological Friedman- 
Robertson Walker expanding-Universe, with weak (but non-zero) particle production. Furthermore 
the D-particle medium can induce an MSW type effect. We have argued previously, in the context 
of bosons, that the so-called flavour vacuum is the appropriate state to be used, at least for low- 
energy excitations, with energies/momenta up to a dynamically determined cutoff scale. Given the 
intriguing mass scale provided by neutrino flavour mass differences from the point of view of dark 
energy, we evaluate the flavour-vacuum expectation value (condensate) of the stress-energy tensor 
of the 1/2-spin fields with mixing in an effective low-energy Quantum Field Theory in this foam- 
induced curved space-time. We demonstrate, at late epochs of the Universe, that the fermionic 
vacuum condensate behaves as a fluid with negative pressure and positive energy; however the 
equation of state has u>f crm ion > — 1/3 and so the contribution of the fermion-fluid flavour vacuum 
alone could not yield accelerating Universes. Such contributions to the vacuum energy should 
be considered as (algebraically) additive to the flavoured boson contributions, evaluated in our 
previous works; this should be considered as natural from (broken) target-space supersymmetry 
that characterises realistic superstring/supermembrane models of space-time foam. The boson fluid 
is also characterised by positive energy and negative pressure, but its equation of state is, for late 
eras, close to tuboson — ► — 1, and hence overall the D-foam universe appears accelerating at late eras. 

I. INTRODUCTION 

During recent years it has been suggested that a certain, mathematically consistent, treatment of flavour mixing 
in Quantum Field Theory could have implications at a cosmological scale Specifically, adopting a Fock-space 
quantization formalism for the "flavour" states [HEl, one can define a new vacuum state which, in the thermodynamic 
limit, is orthogonal to the mass-eigenstate vacuum. This orthogonality, in fact, extends to the entire set of Fock- 
space states constructed out of the flavour vacuum, relative to those constructed out of the mass-eigenstate vacuum. 
It has been claimed [H that the flavour-vacuum formalism, although mathematically consistent, nevertheless leads 
to no physically different predictions from the conventional formalism. However, the authors of || have argued 
that probability conservation is only realised within the flavour state vacuum in quantum field theories with mixing; 
moreover the oscillation probability among flavours is modified, compared to the traditional formalism, by extra terms 
which, although small, nevertheless are in principle experimentally detectable. In this sense they postulated that in 
such cases the flavour vacuum is the physical vacuum. 

It has then been demonstrated [l[ that the vacuum condensate due to fermion particle mixing, evaluated in the 
(physical) flavour vacuum, seems to behave as a source of Dark Energy, in the sense of yielding a non-trivial flavour- 
vacuum energy density. In a series of papers Q , following the initial work of [l[ , it has also been argued that the fluid 
of flavour fermions behaves as an ideal one, in the cosmological sense, with a simple equation of state, which depends 
on the Universe's epoch. 

However, these calculations have been performed in the context of a Minkowski space-time Quantum Field Theory, 
despite the fact that flavour mixing gives rise to a non-trivial space-time cosmological background. A consistent 
treatment, therefore, requires consideration of the Fock-space vacuum in the presence of such cosmological space- 
times, where non-trivial particle production takes place. Moreover, the presence of non-zero vacuum energies, and 
in general of non- vanishing tensor components of the stress tensor of the fermion fluid, indicates breaking of Lorentz 
invariance, except in the case of de-Sitter (or anti-de-Sitter ) vacua, with an equation of state between pressure p and 
energy density p of the form p = —p. For the flavour- vacuum, as shown in Q, and mentioned above, the fermion field 
theory with mixing leads to equations of states that depend on the Universe era, with w — > — 1 approximately only at 
late eras. We should note that, in our opinion, this latter statement has only been argued but not rigorously proved 
in since the relevant calculations have been performed in Minkowski flat space-time. A first such step towards the 
construction of microscopic models that would provide mathematically and physically consistent realizations of the 
flavour vacuum, has been performed in in the context of the so-called D-particle foam model Q, a string/brane 
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FIG. 1: Schematic representation of a D-particle space-time foam model. The figure indicates also the capture/recoil process 
of a string state by a D-particle defect for closed (left picture) and open (right picture) string states, in the presence of D-brane 
world. The presence of a D-brane is essential due to gauge flux conservation, since an isolated D-particle cannot exist. The 
intermediate composite state at t = 0, which has a life time within the stringy uncertainty time interval 5t, of the order of 
the string length, and is described by world-sheet logarithmic conformal field theory, is responsible for the distortion of the 
surrounding space time during the scattering, and subsequently leads to induced metrics depending on both coordinates and 
momenta of the string state. This results on modified dispersion relations for the open string propagation in such a situation 
leading to non-trivial optical properties (refractive index etc.) for this space time. 



inspired model of space-time foam. According to this model, our Universe, after perhaps appropriate compactification, 
is represented as a three brane, propagating in a bulk space-time punctured by DO-brane (D-particle) defects. As the 
D3-brane world moves in the bulk, the D-particles cross it, and for an observer on the D3-brane the situation looks 
like a "space-time foam" with the defects "flashing" on and off ("D-particle foam"). The open strings, with their ends 
attached on the brane, which represent matter in this scenario, can interact with the D-particles on the D3-brane 
universe in a topologically non-trivial manner, involving splitting and capture of the strings by the DO-brane defects, 
and subsequent re-emission of the open string state (see fig. [T]). However, the flavour of the re-emitted state may not 
be the same as that of the incident one, thereby leading to vacuum-induced flavour oscillations and mixing. It should 
be emphasized that, due to electric charge conservation, only electrically neutral matter interacts non-trivially with 
the D-particle "foam" , which is transparent to charged matter 0, |T(| • 

In such a model, the flavour vacuum is regarded as an effective description of the vacuum state for low-energy 
string modes with mixing, as a result of the breaking of Lorentz symmetry locally in space-time, due to D-particle 
recoil during the string-D-particle interactions. Nevertheless, Lorentz symmetry is preserved on the average, in the 
sense that the appropriate vacuum-expectation- values (v.e.v.) of the relevant Lorentz- breaking observables vanish. 
However, this is not the case for the quantum fluctuations of those observables, which may be non-trivial. We consider 
the first-quantised string-theory framework, that describes (perturbatively) the physics of (matter) open strings either 
stretched between the D-particle and the D3-brane world, with both their ends attached to the D3-brane (c.f. fig. [J); 
quantum fluctuations of target-space background fields, in which the string propagates, are induced by appropriate 
summation over world-sheet surfaces with higher topologies (genera). In this sense, a quantum fluctuating D-particle 
in the foam, will be described by such stretched open strings, with at least one of their ends attached to it. 

The structure of the article is as follows: in the next section [H] we review briefly the D-particle foam model and 
discuss the formalism that leads to an induced space-time metric of the form of a conformally flat expanding Universe, 
as a result of the space-time fluctuating D-particle background. In section HTT1 we discuss a gravitational MSW effect 
as a result of the existence of the D-foam, that leads to gravitational-medium-induced flavour mixing and mass 
differences. We also present in that section plausibility arguments for the role of the flavour Fock vacuum and the 
corresponding excitation states as the physical states in the problem. In section HVi in the context of a low-energy 
field theory limit, we discuss the contributions of bosonic degrees of freedom to the dark energy of the brane universe, 
defined as the appropriate vacuum expectation value of the stress-energy tensor between flavour-Fock vacuum states. 
In section [V] we repeat the construction for the fermionic low-energy degrees of freedom, and evaluate the relevant 
contributions to the dark energy of the brane. An important issue arises regarding the choice of the appropriate 
normal ordering that should lead to the physically correct subtraction of the (field-theoretic) ultraviolet divergences, 
in a way consistent with the gravitational MSW effect. This is discussed in detail in section IVTl The reader's attention 
is called at this point to the fact that in string theory there are no actual ultraviolet momentum divergencies. These 
are artifacts of the low-energy local effective field theory which is defined up to energies of the order of the string scale, 
M s , or better the Planck scale (defined as the ratio of M s /g s , with g s the string coupling, assumed weak g s < I). 
Thus, by definition any momentum integral will be cut-off at that scale automatically. The subtraction procedure we 
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are applying to the bosonic or fermionic stress-energy tensors in this field-theoretic context, defines the appropriate 
effective field theory degrees of freedom, accessible to a low-energy observer, and is consistent with the fact that any 
contributions to the vacuum energy should vanish in the absence of D-particle foam effects. In section IVIll we discuss 
the emergence of a momentum cut-off much lower than the Planck scale that arises from statistical arguments related 
to particle production that characterises our expanding background. This is not a sharp cutoff, but rather defines the 
appropriate physical degrees of freedom that lead to significant contributions to the brane- world vacuum energy. We 
also discuss in this section the equation of state of the fermionic vacuum condensate and demonstrate that, for late 
eras of the Universe, it behaves as a fluid with negative pressure and positive energy; however the equation of state 
has Wfermion > — 1/3 and so the contribution of the fermion-fluid flavour vacuum alone could not yield accelerating 
Universes. However, on taking into account the contributions to the vacuum energy coming from flavoured bosons, 
which are natural from the point of view of the (broken) target-space supersymmetry that characterises realistic 
superstring/supermembrane models of space-time foam, and which should be considered as (algebraically) additive to 
the fermion contributions, one may obtain the conditions for late-era acceleration of the Universe. Indeed, the boson 
fluid is also characterised by positive energy and negative pressure, but its equation of state is, for late eras, close to 
Wboson — * — 1 and hence overall the D-foam universe appears accelerating at late eras. Finally, section [Villi contains 
our conclusions and outlook. 



II. D-PARTICLE FOAM INDUCED METRIC 



The target-space quantization of the recoil velocity of a D-particle, Uj, during its interaction with a matter open 
string is achieved [1 11 ] by a genus summation on the world-sheet, which in the case of a bosonic cr-model with a 
D-particle recoil deformation can be cast in a closed form [ll|, [l^] • This yields a stochastic Gaussian distribution of 
the recoil velocities Ui, around a zero average, with a variance a 2 that depends at most on target time, and not on 
the position of the D-particle: 

< m >= , < UiU 3 >= a 2 Sij , a 2 ~ g 2 tl , (1) 

where the time to extends over the capture time of the string by the fluctuating D-particle. As discussed in [T3, E3] , 
for a matter string with total energy po, the capture time is of order t c ~ a'p°. Moreover in [121 ] we also argued, and 
shall review below, that there is a dynamically imposed upper bound scale (cutoff) for the momenta of the particle 
excitations interacting with the D-particles, which is of order of the mass of the particles, essentially (c.f. (|34| . 
below). Even if one considers sneutrinos, their masses (as a result of supersymmetry breaking) may be assumed of a 
few TeV in phenomenologically interesting models, which is still much smaller than a Planck-scale mass of a D-particle. 
Hence the capture time to is small for our flavoured cases examined here. One has to average the relevant expressions 
yielding pressure and energy density of the particle fluid in the flavour-vacuum over such time scales. Because the time 
scales involved are small, such time averages may be replaced, at a good approximation, by the value of the relevant 
quantity to be averaged over a time scale to. This means that on a global scale, on the D3 brane Universe, the quantum 
fluctuating D-particles will fluctuate with an average variance cr 2 (io) 1. On using a dilute gas approximation we 
can take a statistical average over populations of D-particles whose density does depend on the cosmological era of 
the Universe; the D-particle foam could then produce an isotropic and homogeneous (cosmological type) space-time 
background, with non-trivial particle production, in which flavour mixing takes place in a self-consistent way 0], as 
we shall review below. 

For the benefit of the reader, we feel it would be instructive to first review briefly the mathematical formalism 
underlying the quantum-fluctuating D-particle "foamy" space time. The world-sheet boundary operator Vd describing 
the excitations of a moving heavy DO-brane is given in the tree approximation by: 

V D = / {y i d n X i + UiX^nX') = [ Y { (A ) d n X l (2) 

J'dD J'dD 

where dD denotes the boundary of the world-sheet D with the topology of a disk, to lowest order in string-loop 
perturbation theory, m and j/, are the velocity and position of the heavy (non-relativistic) D-particle respectively 
and Yi (A" ) = + UiX°. To describe the capture/recoil we need an operator which has non-zero matrix elements 
between different states of the D-particle and is turned on "abruptly" in target time. One way of doing this is to 
put [IH a O (A' ) , the Heavyside function, in front of Vd which models an impulse whereby the D-particle starts 

moving at A = 0. This impulsive Vd, denoted by V D mp , can thus be represented as 

K np = ^tE / dru t X e (A ) d n X\ (3) 
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where d in the sum denotes the appropriate number of spatial target-space dimensions. For a recoiling D-particle 
confined on a D3 branc, d = 3. 

Since X° is an operator it will be necessary to define (X°j as a regularized operator using the contour integral 



uj — te 



with 



(4) 



where e is a regulator, which, as discussed in [131 ] and will be reviewed below, is linked with a running cutoff scale 
on the world-sheet of the string, on account of the requirement of the closure of the (logarithmic) conformal algebra. 
Hence we can consider 



D £ (x°) = x°e e (x°)=- f 



IS) 



(5) 



The presence of a recoil deformation leads to local distortions in the neighboring space time geometry which can 
be found as follows: let one write the boundary recoil/capture operator V™ p ([3]) in the Dirichlet picture as a total 
derivative over the bulk of the world-sheet by means of the two-dimensional version of Stokes theorem (omitting 
from now on the explicit summation over repeated i-index, which is understood to be over the spatial indices of the 
D3-brane world of fig. [T|) : 
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(6) 



where S £ (X°) is an e-regularised <5-function. This is equivalent to a deformation describing an open string propagating 
in an antisymmetric £^ t „-background corresponding to an external constant in target-space "electric" field, 



Boi 



B; 



, 



(7) 



where the X°S(X°) terms in the argument of the electric field yield vanishing contributions in the large time limit 
e — > 0, and hence are ignored from now on. We remark for completeness at this stage that, upon a T-duality canonical 
transformation of the coordinates [bj ]. the presence of the B-field leads to mixed- type boundary conditions for open 
strings on the boundary &D of world-sheet surfaces with the topology of a disc: 



9^vd n X u + B^ v d T X v 



dv 







(8) 



with B given by 0. Absence of a recoil- velocity H^-field leads to the usual Neumann boundary conditions, while the 
limit where g M „ — > 0, with u\ ^ 0. leads to Dirichlet boundary conditions. 

As discussed in detail in refs. [15|, there is also an induced open-string effective target-space-time metric. To 
find it, one should consider the world-sheet propagator on the disc (X^(z, z)X w (0, 0)), with the boundary conditions 
(j8j) . Upon using a conformal mapping of the disc onto the upper half plane with the real axis (parametrised by r e R) 
as its boundary [l5| . one then obtains: 



(X"(t)X"(0)) = - 
with the non-commutative parameters given by 



clccl.i 



lnr 
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(9) 



\X\t] 



qOI 



lie 1 U 



(10) 



where t is the target time; for definiteness the recoil can be assumed to be along the spatial X 1 direction, i.e. 
^ k\ = k || u\ , ki = &3 = 0. e(r) is the step function having value —1 for r < and 1 for r > 0. The quantity 
Ui = ^ and u c — is the Born-Infeld critical field. Since the space and time co-ordinates are world-sheet fields, 
this commutator is calculated [l5| using the appropriate first-quantised string commutation relations on the world 
sheet. The effective Finsler-type open-string metric [I3|, due to the presence of the recoil- velocity field u (whose 
direction breaks target-space Lorentz invariance) is given by: 



open, electric 
open, electric 



(1 -Ui)T)p V , V,V = 0, 1 

r)^ j v = & H other values . 



(11) 



5 



There is, moreover, a modified effective string coupling [ill, fl^: 

gf = g s (l-u 2 f 2 . (12) 

It should be mentioned here that such metrics have been suggested in the context of a T-dual Neumann picture [Tlj 
of the D-particle recoil process in refs. Q. 

The presence of a critical "electric" field is associated with a singularity of both the effective metric as well as the 
non-commutativity parameter and there is also an effective string coupling, which vanishes in that limit (|12[) . This 
reflects the destabilization of the vacuum when the "electric" field intensity approaches the critical value, which was 
noted in [l8j]. Since in our D-particle foam case, the role of the 'electric' field is played by the recoil velocity of the 
D-particle defect, the critical field corresponds to the relativistic speed of light. This accords with special relativistic 
kinematics, which is respected in string theory, by construction. The critical recoil velocity is of order one, which in 
turns sets the highest order of magnitude of energies of the stretched strings to the mass of the D-particle, M s /g s , 
as announced previously. In this sense, the variances ([T]) are characteristic constants, depending on microscopic 
parameters, such as the mass of the D-particles. 

In the absence of any special knowledge, the most natural choice is to consider isotropic cases of foam, in which the 
induced Finsler metric assumes the conformal form 

open, electric / -i ~2\ „ {-\o\ 

sy£ = f. 1 - j Vnv > ( 13 ) 

for all fi, v = 0, ... 3. This is the case we shall henceforth concentrate on. 

Stochastic quantum fluctuations of the recoil velocity ui, induced by the summation over genera on the world- 
sheet imply, on account of (Q]), a constant metric of the form: 

9^ = (l - er 2 ) Vi* ( 14 ) 

In addition to the quantum fluctuations of the metric as a result of a single D-particle fluctuation considered hitherto, 
one should consider the effects of a statistical population of D-particles, which characterises realistic cases of D- 
particle foam. Denoting such statistical averages over populations of D-particles (as opposed to quantum averages 
over fluctuations of a single quantum D-particle) by <§; . . . one should bear in mind that in general, in situations 
like in fig. [T] the density of D-particlcs in the bulk, which essentially the statistical average depends upon, may vary 
with the cosmological time scale, in the sense that their bulk distribution may not be uniform. Hence, on taking the 
statistical average of the metric (|14p over populations of D-particles, we obtain, in general, a time dependent induced 
metric: 

< 9» v >= gt? = (!- < v 2 > (*)) rinv , (15) 

where <C er 2 ^(t) depends in general on the cosmological time, t, for reasons stated above. In fact, from the conformal 
nature of the induced cosmological metric (|15|) . we see that the time t that appears naturally in our construction is 
the so-called conformal time r\ in standard Cosmology, and indeed the metric effjj* acquires the standard cosmological 
form in the conformal time frame: 

gff = cfa) v > ( 16 ) 

where the scale factor of the D-foam universe is C(rf) = 1— <C a 2 3> (?/). Slightly expanding universes are obtained 
for cases in which the D3 brane moves towards a region in the bulk space (c.f. fig. [J) characterised by a depletion of 
D-particles. Notice that in our construction of small recoil velocities, \ui\ <C 1, for which our perturbative er-model 
treatment of strings suffices, the induced metrics are only slightly deviating from Minkowski space time, and as such 
they constitute good candidates to describe late-eras in our Universe'e expansion. This is the case we shall restrict 
ourselves on in this paper. 



III. D-PARTICLE FOAM AND GRAVITATIONAL MSW EFFECT 



Apart from the induced background space-time P^|) . (fT5|) . the presence of a D-particle foam has other interesting 
consequences for matter flavoured states propagating on the D3 brane world. Specifically, as advocated in [l9j . the 
presence of a fluctuating "medium" of defects in the background space-time, may lead to induced mixing of flavoured 
states, and as a consequence to gravitationally- induced mass differences, analogous to the celebrated MSW effect [20| . 
When neutrinos pass through ordinary matter media such as the Sun, the mass differences and mixing angles acquire 
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parts proportional to the (electronic) density of the medium according to the MSW effect. The difference in the 
(quantum) gravitational case is that the induced mass differences and mixing will now be proportional to the product 
of the density of fluctuating defects in space-time and Newton's constant Gm (which expresses the coupling constant 
of gravity at an effective theory level). In terms of Mp, the four-dimensional Planck mass, we have Gn = j7T (in 

units h = c = 1) . To be precise, it is argued in (l9j . that the gravitationally- induced mass differences among flavour 
states will be of order 

if»L ~ G N (n dcfcct )p (17) 
where (^defect) is an "effective" number density of the space-time defects (D-particles), probed by matter of momentum 

Am 2 

p=\p\. A measure of the weakness of the space-time foam is given by the smallness of the ratio 'I^g"" where m is a 
typical mass scale of the flavoured states. Indeed, for situations in which one has an effective number Af* defects per 
Planck volume, 

AmL m ~W^ M%. (18) 



M P 

The value of Af* takes into account scattering cross-sections of matter with the D-particles and so is much smaller 
than the number of D-particlcs per Planck volume. In the context of the string models we are considering, the four- 
dimensional Planck mass Mp may be different from the string mass scale M s = 1/y/a', which is essentially a free 
parameter to be constrained by phenomenology, according to the modern approach to string theory. 

To ensure Am 2 oam has realistically small mass differences among neutrino flavours, i.e. do not exceed the observed 
values of order at most 1CP 3 eV 2 , one should have sufficiently dilute D-particle foams, such that Af* <C 1. This will 
be assumed throughout this work. A plausible assumption made in [191 ], which, we shall also adopt here, is that 
the induced mass differences are essentially independent of the momentum of the probe, since the effective density of 
defects decreases with the momentum p (the faster the probe, the less time it has to interact with the foam in the MSW 
framework). One should notice, however, that such an assumption is strictly necessary only if the foam-induced mass 
differences are to account for the entire observed mass differences, which are indeed independent of the momentum 
of the neutrinos. This is unnecessary in the case here, since the foam-induced mass differences are only a tiny part of 
the experimentally measured ones, as required by the currently accepted experimental facts [l!|. Nevertheless, due to 
the smallness of Am 2 oam , compared to the typical neutrino mass scales, m, it is convenient to ignore such momentum 
dependence since it suffices to effect the ordcr-of-magnitudc estimates in this article. This will be assumed from now 
on. 

It is important to note that the gravitational MSW effect pertains to flavour mixing induced by the medium and 
does not induce any distortion of space-time per se. According to our discussion in the previous section, it is the recoil 
of the fluctuating defect that induces the background ((15)) . This is an important distinction that should be used later 
on, when we discuss ultraviolet subtractions in our effective low-energy theory of string matter interacting with the 
defects. However, both effects are affected by the density of defects. 

In our model the flavour mixing originates from the fact that, during the capture process of an open string (matter) 
state by the D-particle, the mass m of the re-emitted state might be different from the incident one. In this sense, 
the D-particle "medium" will induce flavour oscillations and mixing, and as a consequence the flavour Fock space 
vacuum is appropriate for quantization of those states, since the "flavour" states are the appropriate physical states in 
this context [7j. Unfortunately at present, our understanding of such a mechanism from a superstring-model, such as 
type IIB string theory, discussed by Li et al. in (To| . is inadequate, due to the non-perturbative nature of the process 
of D-particle-induced mass flips, in a region of strong gravity. Such processes require knowledge of the dynamics of 
D-particles per se, and, unlike the simple recoil/capture processes that do not involve mass changes, cannot be simply 
described by means of (perturbative) world-sheet methods. Nevertheless the time scale involved in such a mass flip 
can be estimated, using stringy uncertainty relations that are independent of the details of the underlying microscopic 
string model. We first note that when a D-particle interacts with a pair of open string states stretched between the 
defect and the D3 brane, representing the capture and splitting process (as in fig. [I}, there is an induced repulsive 
short-range potential V, calculated by means of appropriate world-sheet annulus graphs in [2l|, following techniques 
developed in [13] • The relevant parts of the potential for our discussion in this paper are of the form 

7r a' u 2 . . 

12 r 6 

where u is the relative four-velocity (i.e. u 2 — jzrp in units of the speed of light in (Minkowski) vacuo c — 1 and v 
is the 3-velocity) between the D-particle and the D3-brane and r is the distance between the defect and the brane 
along the transverse directions to the brane world (c.f. fig. [TJ there is no potential for D-particle motion parallel 
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to the brane (22J) 1 . In our case, with a fundamental string stretched between the two, the order of this distance 
can vary typically from that of the string length V" = 1/M S to a (much smaller) characteristic minimum one 
L m in ~ tHo/ *C v / ai' where to is a typical mass scale of the stretched light string state representing the flavoured 

states HH, HI]- The corresponding magnitudes |V| are |V| ~ M s u 2 and |V| ~ (y=-) M s u 2 . In this work we will 
restrict our attention to the simple case of two dominant mass eigenstates, with masses rrii,i — 1,2; one may then 
take to ~ \{mi + m-x). It will be convenient to collectively represent the effects of these two extreme cases for the 
potential in a single formula: 

|V|~ (^-X M s u 2 , (20) 

with q = 0(3) for the former (latter) case. 

A D-particle induced mass flip of a flavoured matter string state, such as a neutrino, will be between masses 
separated by terms of order (fTT|) . (fTS)) in our scenario. Since the recoil contribution of the D-particle cancels out 
on average, one has correspondingly a momentum conservation for the neutrino state interacting with the defect. In 
the non-relativistic limit of the D-particle recoil velocities |u| ~ \v\ -C 1 energy conservation and in a first quantized 
weakly coupled g s < 1 string theory framework, energy conservation implies: 

M D + < (p 2 + To! 2 ) 1 / 2 >=< (p 2 + TO2 2 ) 1 / 2 > +M D + - M D < u 2 > +C(u 4 ) , (21) 

where Mp> = g ~ Mp is the D-particle mass (assumed to be of order of the four-dimensional Planck mass Mp). 
This allows us to estimate the recoil velocity fluctuations during a mass-flip process. The above relation expresses 
an average over both quantum fluctuations of the recoil velocity at the individual D-particle level and over D-particle 
foam populations. This averaging is denoted by <C ... 3>. On denoting the neutrino energy difference by AE, one 
may estimate from (|2"Tj) : 

M s < u 2 >~ 2g s AE . (22) 
Hence, the recoil velocity fluctuation accompanying a mass-flip induces in turn a fluctuation in the potential (|19|) of 

order (for r ~ (fny/a'^j y/cV, with q' = or 1, corresponding to the two characteristic scales of the intermediate 
string states discussed above) : 

AV~2 0^ g s AE. (23) 

In the laboratory frame from the saturation of the energy-time uncertainty relation an energy fluctuation will 
imply a life-time ^mass-flip for the intermediate string state associated with the mass flip. It suffices to consider the 
low-energy quantum mechanical version of this uncertainty, which yields (in units h = 1) 

1 / to \ 9 1 

A< — -*p~AV~UJ 2g-AE> 9 = ' 3 ' (24) 

We will denote by A£ cap t U re the time during which an intermediate string state, stretched between the D-particle 
and the D3-brane world, grows from zero size to its maximal one permitted by the stringy time-space uncertainty 
relations, and back to zero size. To have the possibility of mass-flip in string theory, At^ lass _ flip must be longer than 
the time involved in capture, At capt uro- As discussed in [Io| . the capture time is of order 

a a 'P° 

^^capturc ^° Z ^ , V / 

1 — U z 



1 In the supersymmetric model of [2l| . the brane is a D8-brane, from which a D3-brane can be obtained by appropriate compactification. 
The form of the (repulsive) short-range interaction, and in particular the dependence on u and r, are insensitive to such compactification 
details, but do depend on the dimensionality of the interacting branes [2^1 . On the other hand, the precise form of the numerical coefficient 
depends on the details of the construction and on the presence of other branes and orientifold planes (the latter ensuring dynamical 
compactification of the bulk space in the model). However, for our order of magnitude estimates such issues are not relevant. Moreover, 
the presence of orientifold planes in the construction of \'2\\ can cancel velocity independent terms in the potentials. Such terms are also 
irrelevant since we are only interested in potential fluctuations induced by the velocity fluctuations, see eq. (1231 below. 
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where u 2 is the D-particle recoil velocity, and p° ~ (p 2 + to 2 ) 2 is the energy of the incident string state. 

The delays ([25]) are consistent with the time-space uncertainty relation At AX ~ a' , characteristic of string the- 
ory [23l ] and actually saturate it. They can be computed rigorously within superstring theory in D-particle backgrounds 
by evaluating the relevant scattering amplitudes and looking at backward scattering contributions flOl . flq . In what 
follows we shall only be interested in contributions to leading order in small quantities, and so, when using ([25]) . the 
recoil velocity u will be ignored. The capture time (|2"5"|) does not involve mass flip, and it is associated simply with 
capture and re-emission of an open string state by the D-particle. 

We reiterate that, in order for the mass-flip process to be feasible within a string theory model, it is necessary that 
the capture time (|25|) . which saturates the stringy time-spaced uncertainties, is shorter than, or at most equal to, 
the mass- flip time (|24p . Otherwise, mass flip does not take place. Such a requirement, then, implies an upper bound 
for the spatial momenta of the states that can possibly undergo mass flip, i.e. for the momenta associated with the 
Fock-type flavour vacuum 0, 0] 

First let us consider the case q — 0. It is easy to see from ([24]) and ([25]) that for all momenta that define an 
effective-low energy theory, i.e. momenta smaller than the Planck scale p < M s /g s , the condition 

AC-flip > ^capture (26) 

is comfortably satisfied, since to violate it requires energy scales 

g s 2AE g s 

For a ratio £ of typical neutrino mass scale to to momenta, with £ either small or large, compared to unity, this 
inequality can only be satisfied for mass differences 5m 2 = \m 2 — m|| > M 2 / g s . This requirement is physically absurd, 
and thus yields no physically sensible constraint on the neutrino mass differences for the case q = 0. 

However, on taking into account the case q = 3, the condition (|26[) is satisfied for all momenta smaller than the 
cutoff-scale M s /g s provided 



Sm 2 < m 2 (28) 

which is a strong constraint for the foam-induced mass differences. 

The bound ([28]) leads to very small foam-induced mass differences for large string mass scales, M s , close to the 
Planck scale Mp, while one can get mass differences of the order of the observed ones for M s of the order of TeV, and 
g s <~ 10~ 16 , such that M s /g s — Mp ~ 10 19 GeV. Compactification details, of course, in phcnomcnologically realistic 
string/brane models may affect such estimates seriously. The reader should notice that in the bound ([28]) there is a 
theoretical uncertainty at most of order O(10), due to the uncertainty in the numerical coefficients in the potential 
(|19p . which depend on the details of the microscopic model, as already mentioned. 

Hence mass flip processes, at least from the point of view of stringy uncertainties, are consistent physical processes 
that can take place for neutrino energy differences of physical relevance. A microscopic understanding of these 
processes in detailed realistic superstring/supermembrane models is, of course, still pending, and hence we can only 
give here plausibility arguments on the existence of such processes. 

It should also be noticed that, for small masses over momenta, and on assuming, for concreteness, only a D-foam- 
induced mass difference ([T8]) among flavours, one may estimate from ([2T]) the average stochastic fluctuations of the 
D-particle foam recoil velocities: 

« u 2 » = « a 2 »~ g, A "f oam ~ W* <C 1 . (29) 

M s p 

Here we should remember that M* denotes the effective number of D-particle defects contained in a Planck volume. 
The reader should then notice the form similarity of ([29]) with the gravitational MSW-like relation (]17p conjectured 
in ref. [l||. The order of the estimate, of course, may change significantly if the induced mass differences among 
neutrino flavours due to the foam constitute only a small percentage of the physically observed one ( as is most likely 
the case) even if the D-particle foam is physically relevant [l9| . 



IV. LOW-ENERGY BOSONIC FIELD-THEORY MIXING AND FLAVOUR VACUA 



In this article, we shall discuss mixing of field theory excitations, induced on both bosonic Q and fermionic 
excitations of strings by D-particles, during their topologically non-trivial interactions (spiltting/capture/re-emission) 
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with strings. However, the mixing will be discussed in the presence of a slightly expanding universe (|15[) . induced 
on global scales by time varying populations of quantum fluctuating D-particles, as discussed above. The presence 
of both fermionic and bosonic "flavoured" field theory excitations of strings, finds a natural application in the case 
of superstrings and super-D-branes, which is the ultimate physical theory we have to consider. Indeed, even if 
supersymmetry is eventually broken in target space, the partners (e.g. sneutrinos) do exists, and in the case of flavour, 
their mass differences might be of the same order as the original particles (corresponding to flavoured neutrinos in our 
example), despite the fact that the mass differences between particles and their supersymmetry partners might be at 
least a few TeV, due to the broken supersymmetry. The only caveat in our mathematical construction is that it is 
based on bosonic string theory, where the resummation of leading modular divergencies in the case of recoil- velocity 
deformed cr-models, is possible Unfortunately in the supersymmetric case, which would necessitate world-sheet 
supersymmetry as well, such a resummation is not possible at present (24[. Thus for the (realistic) case of (broken) 
supersymmetric D-foam we could only assume that the conclusions drawn from the bosonic case, regarding stochastic 
fluctuations properties of the foam, are sufficiently robust to be extendable here. 

We commence our discussion with a review of the bosonic case. The (l+l)-dimensional bosonic case has been 
discussed in detail in 0] and will not be repeated here, apart from pertinent information needed for completeness of our 
discussion. It has been shown there that the vacuum condensate in the case of scalar fields, taken as a representative 
example, behaves as a fluid with w ~ —1 once the MSW effect is taken into account. The mixing/expansion is 
considered in a (1 + 1) dimensional framework, since the recoil of the D-particle has been taken parallel to the motion 
of the bosonic excitations, assumed along one spatial direction, say X 1 (c.f. (|f I[) . In view of our isotropic foam 
situation, considered here, (|13[) . (TT5)) our findings should carry forward to the full (3 + I ) dimensional case. 

For a scalar field <f>, the stress-energy tensor is: 

T*° s [4>] = \ (D^D^ + D^D^) - g^L 

The 1 + 1 dimensional version (with g^ v — C^Tj^v in the conformal frame (c.f. (jllj) ) was considered in 
was shown that the only non-trivial components are the diagonal ones: 

T*° s [<t>] = (d^f + id^f-Ceff^m 2 ^ 2 (31) 

The reader should have noticed that we used the symbol C e // in the mass term for the scalar field and not simply 
Cm 2 4> 2 . The quantity C e // contains both the effects of the expansion and the MSW effect and plays the role of an 
effective scale factor; this is in a similar spirit to standard effective field theories of inflation [25j | with interactions 
of massive dark matter particles to the inflaton field. The dispersion relations of such massive particles, and the 
associated stress tensor components, include the effective scale factors, as above. We shall come back to this point 
with more details when we discuss the fermionic case, where an entirely analogous situation applies. 
The expressions (|3Tj) readily generalize to 3 + 1 dimensions as follows: 

3 

4o s M = (M 2 +E^) 2 +Ceff(v)m 2 <f? (32) 
i=i 

3 

2*»M = (d v <f>) 2 + ]T(cM) 2 - C eff (r,)m 2 <f> 2 + 2{d x ^f 
i=i 

where the sums over spatial three-dimensional indices are explicitly denoted for clarity. Since the form of this term 
does not change when we go from the 1 + 1 dimensional analysis to the 3 + 1 one, the conclusion that w ~ — 1 is valid. 
As discussed in Q , the appropriate normal ordering (subtraction) in our case of D-particle foam has to remove any 
terms that do not vanish in the limit where the variance of the D-particle fluctuations vanishes, <C a 2 ^> — ► 0. 
We further assume that at late eras of the universe the D-particle fluctuations are weak, so only leading order terms 
in an expansion in powers of <C <r 2 3> are kept. 

Hence, after the appropriate subtraction, discussed further in [7|, one arrives at 

,<0 |: T^ s [4>] :\ 0) f = f (0 \: C e// (r/)m 2 ^ 2 :| 0)/ 

,<0 |: T*r[4>] :\ 0) f = - f (0 |: C eff ( V )m 2 ^ 2 :| Q) f . (33) 

Any extra contributions acquired due to the higher ((3 + 1)) dimensionality, as compared to the (1 + l)-dimensional 
case of ref. 0, are common to the two components; therefore the equation of state w w — 1 holds also in the 3 + 1 
dimensional bosonic field case. 
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To recapitulate, for weak D-particle foam, the energy density of the bosonic fluid in the flavour-vacuum formalism 
is positive, while the pressure is negative, and the equation of state is consistent with a cosmological constant. There 
is however an important issue to be addressed here. The expressions in (f33|) . which involve integration over momenta, 
are formally ultraviolet divergent 0]. Hence in the low-energy field theory limit, the momentum integrals need a 
momentum scale cut-off k max . It should be emphasised that in string theory there are no actual ultraviolet momentum 
divergencies. As already mentioned in the Introduction, these are artifacts of the low-energy local effective field theory 
which is defined up to energies of the order of the the Planck scale M s /g s . Thus, by definition any momentum integral 
will be cut-off at that scale automatically. The appropriate effective field theory degrees of freedom, accessible to 
a low-energy observer, are defined by the subtraction procedure that we are applying to the bosonic or fermionic 
stress-energy tensors in this field-theoretic context. Moreover any such contributions to the vacuum energy should 
vanish in the absence of D-particle foam effects. It is in this sense that a cut-off k max is used in the model. 

In 0], such a cutoff scale, which is, however, much smaller than the Planck mass, was determined dynamically, by 
considering particle production. The result of the (1 + l)-dimensional case of Q has yielded the cutoff scale 

*W ~ , (34 ) 

where m^, i = 1,2 are bosonic eigenstate masses. This is a result of the fact that the particle production falls off 
with the momentum, in such a way that the vacuum is populated significantly by flavoured bosons for momentum 
scales below (|34[) . and thus it is in such regimes of four- momenta that the condensate becomes significant. For small 
mass differences, compared to masses, which we assume throughout our works 2 , we may use the parametrization [7]]: 

(±) 8m + _ 

TTV ' — 777 ± — — , 777 = 777l , "7 = 7772 

. _ , i (M 2 

"max ^ ?77 i n , 

8 771 

Sm = 777l — 1T12 <C 777 . (35) 

Upon inserting the cutoff in the momentum integrals and performing the appropriate subtractions (|33[) . based on the 
metric (|15p . we arrive at the estimate for the boson- induced flavour- vacuum energy density 

Pbosons ~ sin 2 6» < cr 2 (t ) > (W) 2 (36) 

in the case of predominant two-flavour mixing, which we restrict ourselves here for concreteness and brevity. The 
estimate requires time averaging of Q33[) over small capture times, to, which in our model has been estimated to be 
close to zero Q. This fact makes oscillatory terms, that may appear in the components of the stress tensor, negligible. 
Notice in our case the extra suppression factor <C c 2 3> (to), which for the cases of weak gravitational foam is smaller 
than one, as compared to the flat-space time flavour vacuum case of |l|, |6|. In view of (|36[) bounds on <C u 2 3> may 
then be imposed by cosmological considerations, given the order of magnitude of the dark energy at present eras of 
the Universe, observed today. It is important to note that, in our approach, we assume that the relative motion of 
D3-brane worlds in the bulk at present eras is such that the associated supersymmetry breaking due to brane motion 
and the pertinent contributions to vacuum energy are negligible compared to the flavour- vacuum ones (|36p . This is 
an assumption that holds also for the fermionic contributions. 



V. LOW-ENERGY FIELD-THEORY MIXING AND FERMIONIC FLAVOUR- VACUA IN 

(3+l)-DIMENSIONS 

As we shall discuss below, as far as the equation of state is concerned, a crucial difference appears when one 
considers fermionic low-energy field theory excitations in the context of flavour vacua. When the normal ordering 
procedure is applied, the resulting equation of state is of the form > w > —1/3, a range insensitive to the cut-off. 
Thus, fermions alone, cannot lead to a present-epoch acceleration of the Universe through this mechanism. However, 



2 We remark that, even if the bosons refer to sneutrinos, which have heavy masses due to target-space supersymmetry breaking, the 
relative mass differences between mass eigenstates may be assumed sufficiently small, since the mass differences are independent on 
supersymmetry, especially if, according to our D-particle foam model, these mass differences originate from foamy interactions, and 
hence are quantum-gravitational in origin. 
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in our D-particle supersymmetric foam, where, for reasons stated, one has both bosons and fermions as a result of 
broken target-space supersymmetry, the contributions to the equation of state of the flavour vacuum from flavoured 
bosons (e.g. sneutrinos) can lead to a current era acceleration, by affecting the equation of state appropriately. 

Before presenting the details of our calculation, let us briefly summarize the formalism for the mixing of two 
fermionic (spin 1/2) flavours in a Minkowski space-time background: two flavoured fermions ipe(x) and tp^ix) can be 
constructed from two free Dirac fields ipi (x) and ipzix) with definite masses mi ^ m2, by means of the relation 

ip e (x) = tpi(x) cos 9 + ip 2 (x) sin 9 , . 

ip^{x) = —ipi(x)sm.Q + il)2{x)cos9 ^ ' 

It has been shown 0] that in quantum field theory (QFT) it is possible, in a finite volume, to define rigorously a 
unitary operator that behaves as the generator of the mixing transformation for fields: 



4> e (x) = G\(t)Mx)G e {t) 
M x ) = Gl{t)Mx)G s (t) 



(38) 



This can be written as 



G e (t) = cxp 



(39) 



9Jdx ($(s)&(a;)-$&(aO) 
thus allowing the definition of a flavour vacuum as the state 

|0>/ = £t(t)|0) (40) 

with | 0) the mass eigenstate vacuum used in the quantization of the field theory of ipi(x) and xf>2(x). We should 
recall that in the thermodynamic, infinite-volume limit, as in the bosonic case, the vacuum states | 0) and | 0) f are 
orthogonal. Indeed this is the case for all the Fock-space excited states constructed out of these vacua However, 
as discussed in 0] and also below, in our string-inspired effective theory, flavour boson or fermion states, constructed 
out of | 0) f, will exist up to a dynamically determined momentum cutoff scale. Above that scale, the pertinent states 
will be constructed out of the mass-eigenstate vacuum | 0) . This is one way that our microscopic approach differs 
from that of H H. 

We wish now to study how the flavour vacuum expectation value of the stress-energy tensor varies in a weakly 
curved space-time background. For our model of D-particle foam, this is crucial, since the mixing phenomenon and a 
non-trivial, Robertson- Walker type, space-time background (|15[) are induced as a result of the quantum fluctuations 
of the D-particles. Thus, self consistency requires an analysis of the flavour vacuum in the presence of a non-trivial 
curved space-time, where particle production takes place. To this end, it is of interest to first consider a free theory 
of two fields of definite mass in a curved space-time, that will be regarded as a classical background. Then, under the 
assumption 0] that the metric is asymptotically flat at early times, one can introduce the mixing (and therefore the 
flavour vacuum) at t — > — 00 in the way presented in the literature (2j. In our model of D-particle foam (c.f. fig. [l}, 
the assumption of asymptotically flat space times at early times, t — ► —00, can be justified in a scenario in which the 
D3-brane worlds at t — > —00 finds itself in a bulk region depleted of D-particles. As the cosmological time elapses, 
the D3-brane world may move into a densely populated bulk regions of D-particles, and subsequently exit them at 
asymptotically long times in the future, t — ► 00, in such a way that an interpolating cosmological space-time between 
— 00 and +00 may be induced as a result of such D-particle configurations. For current eras of the universe, we may 
assume that the density of D-particles is low enough so that there are only slight deviations from the Minkowski 
vacuum. 

Adopting the Heisenberg picture [2j, |7| , one can study the vacuum condensate induced by the mixing. Hence the 
flavour vacuum expectation value (vev) of the stress-energy tensor operator is evaluated for a time-independent flavour 
vacuum state and a stress-energy tensor that evolves with time. Therefore, we are going to construct the stress-energy 
tensor operator for a theory with two free fermions in a curved background and evaluate its vev, considering the flavour 
vacuum defined in (|40[) and (|39[) as the vacuum, for t — > —00. 

As has been well explained in the literature [26[ , the classical theory for fields with spinorial structure is generalized 
in curved space-time through the vierbein formalism. The expression for the stress-energy tensor for a free 1/2-spin 
field reads 

Tp V (ip) = -g^C + i ($7( /1 D 1 ,)V' - D^Tp^^ip) (41) 

where C denotes the Lagrangian of our theory, 7^ are the generalized 7- matrices defined by {7 M ,7i/} = 2<7 Mt/ , D M is 
the gravitational covariant derivative and ip = iipj a , and 7 is the temporal component of the ordinary 7-matrices 
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in the tangent plane (denned by {j a ,j b } = 2i] ab , a,b — 0, . . . 3 tangent plane indices in the vierbein formalism). 
The stress-energy tensor for two free fields is simply T tlv {^i,^ 2 ) — T^fyi) + T^^)) therefore we can focus just 
on Tfiu^tp). Assuming that the flavour mixing does not affect the homogeneity and the isotropy of the universe we 
can consider a Friedmann-Robertson- Walker (FRW) metric in conformal coordinates {rj, x}, = C^rj^, with 
C(rj) > 0. In this case, adopting the convention rj^ = diag{—l, 1, 1, 1}, the stress-energy tensor becomes 



Tp V (i>) = -Cn^C + ipl—j^d^ + — -^7 (|1 [7o,7 l ,)] V + h-c. (42) 




An important remark on notation should be made at this point: from now on (i.e. in all subsequent formulae) 
the contraction of the space-time indices fi, v . . . is understood with respect to the Minkowski part of the metric, 
i.e. A^B 11 = rj^A^B" ', since the scale factor C(r)) has been factored out appropriately. The notation regarding the 
contraction of tangent-space indices a, b . . ., remains as before, the contraction involving the tangent-space Minkowski 
metric rf h . 

The quantization of this theory has been specifically carried out by [13] and discussed more generally by 28]. 
According to [27| . the quantized spinor field can be written as 

^,x) = (— =) ^^''(Ml^'^M)^^^'^'') (43) 

a, b = ±1 

where L is the parameter that enters the periodic boundary condition ip(rj, x + nL) = ip(r), x) (n being a vector with 
integer Cartesian components), 07(77) = \/p 2 + m 2 C{r\), v^ a,b \p,r]) is a spinor defined by 

( p, 7?) EE {p/ \/CM) 

(~iay/p 2 + m 2 7° + iaj ■ p + m)v {a - b) {p) = (44) 
« (o ' b)+ (^ (a '' b 'H$ = <W V 

and a,( a > b \p, rf) are operators such that {a^ a - b \p, 77), a^ a ,b ^(q, 77)} = S a , a '^b,b'Sp y q. In order to introduce the mixing at 
early times, we first have to define our Fock space for 7/ — > —00. The relevant operators are A^ a ' b '(p) = a( a > b '(p, —00) 
with {A( a > b \p), A^ a ' b ^(q)} = Sa,a>5b.b'5p,q- Starting with these operators, we can define a Fock space for 77 — > —00 
following the usual prescriptions of minkowskian QFT. Moreover, if we assume that for our conformal scale factor 
behaves as C(— 00) — > 1, we can introduce the mechanism of the mixing in the way that has been explained above. 
Therefore the flavour vacuum will be defined by 

I 0)/ ee G+(-co) I 0) (45) 

keeping in mind that, since the space-time is asymptotically flat at early times, the fields ip behave in the usual 
relativistic way: 

7^(77, 2) = iV A {a ' b) {p)v {a ' b) (fi e ™(p- s -Vp 2 + m2 v) (46) 

P 

a. b = ±1 

for 77 — > —00. (N.B. the mass-eigenstate index i is being omitted here for brevity.) We want now to evaluate the 
evolution in (conformal) time of the flavour vacuum expectation value of the stress-energy tensor. Noticing that the 
flavour vacuum has been defined in (|45[) by means of the operators A^ a,b \p), whereas the fields ipi are defined at any 
time in terms of the operators d^' 6 - 1 (p, 77), we have to find explicit relations between these two sets of operators that 
hold at any time 77. This can be achieved by noticing |27j that 

a (a - b) (P, V)=J2 D c v)A {c ' abc) (acp) (47) 



a— ±c 



with D% (p, rf) defined through the equation 



D$(p,v) = C +a'£drf\ ^ ™ )2 e^I^'D^ a , )M (48) 
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with a, a' = —1,1. Hence, at any time 77, the solution of the equation of motion (|43[) can be written as 

#(r,,%=(—±=^J ^i^'^^^fe^^^^^^e^-^^^^ (49) 

a, fa, c = ±1 

The functions D£(p, 77) encode the features of curved space-time, but an explicit solution of equation (|48|) for a generic 
conformal scale factor C(rf) is not known. In the rest of the work we will keep them in their implicit form, with the 
understanding that explicit formulae have to be evaluated case by case. In particular, in our case of D-particle foam, 
in which for late eras of the Universe the conformal scale factor of the induced metric (|15p is close to one, for small 
variances <C er 2 3> <c 1 of the D-particles, we can evaluate approximately such functions, in an appropriate expansion 
in powers of <C o 2 

These tools are sufficient to show that, in the continuous limit (L — > 00), we have 



f (0 I Tu^ufa) I 0)/ = {0| T«(V>i,Va) I 0) +sin 2 0/ dp [V 2 (p) {T u ( V ,p, mi) + T u (v,p, m 2 ))] +0(sm 3 9) (50) 

Jo 

and 

f (0 I T 00 (^i,V2) I 0)/ = (0 I Too^i,^) I 0) + sm 2 9 dp [V 2 (p) (T 00 ( V:P , mi) + T 00 (r),P, ma))] + 0(sin 3 0) (51) 

Jo 

with all the other components vanishing. Above we used the notation 



Tii{r),p,m) = --TT-TTr 73— (1 - Pi | — |-D_i| ) H #e 



3(2tt) 2 w(r?) \ v 111 p 



Too^p.m) = ^^(^ycM^-^r 1 ! 2 -!^-!! 2 ) 



(52) 



and 



_ , _ ^/p 2 + m 2 -y/p 2 + m 2 — p 2 — mim 2 



^(p) - v " 7- V " 7. = g — ■ (53) 
y p + m\ y p + 

As in the boson case |Y| , we have also introduced a cutoff fc max in the four-momenta, since all expressions are affected 
by ultraviolet divergences. 

Before discussing these divergences let us stress that in case of 9 = and/or mi = (i.e. in the absence of 
mixing), the function V 2 (p) vanishes and therefore /(0 | T^fyi,^) \ 0) f — (0 \ T tiu (ipi,ip 2 ) | 0), consistently with 
what we would expect if we set tpi = -02 in (|4"0"]) and (|39p 0. It is important to remark once more that, in the context 
of the Heisenberg picture we have adopted here, the flavour vacuum condensate V(p) 2 , which encodes the structure 
of the flavour vacuum in terms of Fock states in the mass representation [2] , has been computed in the asymptotic 
past (77 — > —00), when the induced space-time (|15j) is assumed flat (assuming a uniform D-particle background at 
77 — > — 00, with no significant curvature effects. It is understood that such assumptions are highly model dependent, 
and hence the considerations here are specific to the initial conditions). In our Heisenberg picture, the term V(p) 2 
does not evolve with time, and as such, its contribution will not be involved in any normal ordering. 



VI. NORMAL ORDERING 



Let us now focus on the ultraviolet divergences. These infinities, in local field theories, are conventionally removed by 
renormalization. In flat space-time this has been achieved by a suitable normal ordering. In the case of a conventional 
local field theory in a time dependent metric background, the renormalization 28] would involve state dependent 
counter terms which, in a covariant procedure, can be tensorially constructed from the metric tensor. As already 
mentioned, in the present field theory model, which is the low-energy limit of a string theory involving D-particle 
capture of stringy matter [3] (c.f. fig. [T]), the procedure of normal ordering is dictated by the underlying microscopic 
physics; the subtraction procedure, in particular, has to be such that in the limit of the absence of D-particles and 
their fluctuations, <C o 2 3>— > the mixing phenomenon should disappear. Moreover in our D-particle foam model 
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one needs to distinguish two effects, as far as the structure of the underlying space-time is concerned. The first effect 
concerns a background space-time, over which propagation of low-energy matter excitations (fermions or bosons) 
takes place. The background space-time in our case of D-particle foam has been argued to be obtained from quantum 
fluctuations of individual D-particles, which in the case of first quantised string framework are due to a summation 
over world-sheet topologies [ll| , upon (statistically) averaging over populations of D-particlcs on the D3-brane world 
(c.f. fig.Q}. This leads to a background space-time metric of the form (|15[) . 

The individual MSW interactions of the flavoured matter excitations with the D-foam background [l9| . produce 
extra back-reaction local fluctuations on the space-time structure. They do not cause metric distortions, as already 
mentioned, but affect the particle mode's energy-momentum dispersion relations. This parallels and is in a similar 
spirit to the standard result when one considers particle production at the end of inflation [25[ . Hence in the effective 
QFT we have to take into account interactions of the neutrinos and the D-particle medium. In this respect, the scale 
factor C(rj) appearing in the above formulae for the fermions, should be considered as representing the background 
space-time. The energies uj on the other hand will contain an 11 effective" scale factor 



In the dispersion relation 



c',.,r(//i=t'(//) + Ac'(//l^t'(„)(l+^M). (54) 



^cff = Vp 2 + m 2 C cS (ri) - VP 2 + m*C(r,) + — ^ ^£^= , (55) 



2 \Jp 2 + m 2 C(r]) 



to leading order in the approximation |AC| <C 1. The AC is the MSW contribution. For the MSW scenario (c.f. (|18[) . 
(|2"9")0 . one has the estimate 

m 2 AC(r 1 ) ~ C(ry)Am? oam ~ C (n) TV* M P p ~ C(n) « a 2 » g^ 1 M s p . (56) 

Such dispersion relations, that take proper account of the non-trivial interactions of the matter probes with the 
D-particles in the foam, should replace the free-particle dispersion relations considered in [27[ and used so far. 

Notice that in our model, in general AC would depend on both momenta and position of the mode. For the purposes 
of this work, we are only interested in an order of magnitude estimate of the induced dark energy contributions. We 
do not need to specify the precise form of the induced distortion AC, apart from the fact that in the context of our 
model we know that this will be proportional to the stochastic fluctuations of the recoil velocity, o- 2 (rj), at the time 
of the interaction. For sufficiently small <C <y 2 S> one may ignore the momentum dependence of AC (n) in (156p . On 
this basis, one can replace u> in Eqs. (I52|) by ui c g (|55p . leaving the initial-time condensate V(p) 2 intact. In a similar 
vein, eq. ([48]) . becomes: 

d&m = e+'£*0AWs^^.H^i>w (»,/)) (57) 



with a, a' = —1,1 and uo(ri) = \J p 2 + m 2 s (rf)C{rj) : and m 2 s C(rj) = rn 2 C c s(jf). 

After these necessary preliminaries, we are now in a position to discuss our subtraction (normal ordering) procedure. 
The latter, will be defined in such a way that in the absence of any MSW interaction of the fermion matter with 
D-particles, the stress tensor should vanish. Notice that the MSW interactions, which are proportional to the density 
of defects, are disentangled from the space-time background C(rf), in the sense that they contribute to mass shifts 
but they do not induce space-time distortions. However, in view of (I29p . in our picture, the induced mass shifts are 
also proportional to the stochastic fluctuations of the recoil velocity of the D-particles, which affects the space time 
background, c.f. (fT5|) . Hence, by subtracting the MSW-like terms from the stress tensor, no further subtraction would 
be necessary in order to ensure that in a Minkowski space time the stress tensor would vanish. Replacing, then, u 
by u e g in Eqs. (|52|) one should use (|55|) to leading order in a small- AC (77) expansion to determine the energy density 
and pressure of the fermion fluid. 

Since in our D-particle case, the change AC is proportional to the variance cr 2 (r?) <C 1 of the recoil- velocity fluctua- 
tions in our weak space-time foam background (c.f. (|29p ). only leading order contributions proportional to a 2 should 
be taken into account. We observe that the terms involving operators become irrelevant to this leading order, 
and so the latter should be replaced by their flat-space-time counterpart, which vanish (c.f. (|48|) . ([57]) ). As a result, 
Eqs. (|52|) become: 



1 



2 

■E 



3(2tt) 2 t!(p 2 + m 2 C{n)) 
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2 



T 00 (r],p,m) :« — -— 5 — AC(t7) 



2 



( 2 ^) 2 £i (p 2 + mfC(r;))^ 

(58) 

where the approximate sign indicates leading orders in AC. Since, to leading order in a 2 , it is expected, rather 
generically, that <C AC (7/) ;g>cx<§; er 2 (r;) the background scale factors C(rj) in the above relation can be replaced 
by constants (i.e. flat space-time). Moreover, upon taking the statistical average <C . . . 3> of (fB"8"|) over D-particle 
populations, to a good approximation for the weak space-time foam situations of interest, any momentum dependence 
of <C AC (77) ^> disappears; hence the latter quantity can be taken out of the momentum integrals in ([50|) . (|5"Tj) . This 
will be understood in what follows. 

On recalling that for a relativistic fluid Too/C(rj) represents the energy density, and Tu/C(rf) the pressure, we can 
easily see from (|58p . (|50[) . (1511) that our fermionic vacuum condensate behaves as a fluid with negative pressure, and 
positive energy density, with an equation of state that satisfies —1/3 < w < 0. This is the result of the opposite 
powers of ui (and hence oj e g, according to our discussion above) appearing in a specific way in the pressure and energy 
expressions dSOj) , dSTJ , (f52|) . 



VII. EQUATION OF STATE, VACUUM ENERGY ESTIMATES AND DYNAMICAL MOMENTUM 

CUTOFF 



To determine the precise value of both the energy density and pressure, and hence the equation of state, it is 
necessary to have knowledge on the momentum ultraviolet cutoff fc max used to regulate ultraviolet infinities in flat 
space-time. To leading order in the small expansion AC (77), such flat space-time approximation for the evaluation of 
the cutoff function proves sufficient. As discussed in Q, a dynamical cutoff function appears if one considers particle 
production due to the flavour vacuum. In flat space-times, the particle number is given by: 



,(0 I N e (p) \0) f = sin 2 9 H(p) (1 - cos [fa + lo 2 ) rj\) + 0(sm d 9) + 0(a 2 ) 



with 



S(P) 



((uj 2 -m 2 )(wi +m 1 ) -p 2 Y 
2luilu 2 (u2 — m 2 )(wi + mi) 



(59) 



(60) 



and u>i = \/p 2 + mf, i = 1,2. The p dependence is essentially determined by the behaviour of 5(p). For large p 




FIG. 2: The function H(p) is plotted in the range p £ [0, 10] for the values mi = 1 and 7712 = 2 (in arbitrary units, just for 
illustration purposes). Notice that there is a maximum Q], corresponding to the point y/m\mi = v2 3 in our arbitrary units. 



(compared to masses) the function falls off with an inverse forth power of momentum (for a plot of this function vs. 
momentum, p, see fig. [5]): 



S(P) 



1 (m\ — m 2 ) 2 (777,1 — to 2 ) ((toi + m 2 ) + 2(m\ + ra 2 .)) 



2 p 2 8p 4 
and so there is a scale (determined by the ratio of the two terms in the above) 

fc ~ - \J 3m 2 + 2m 1 m 2 + 3m 2 , 



(61) 



(62) 
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which is a plausible cut-off scale in momenta p. Thus although there is no sharp cutoff, nevertheless, the flavour 
vacuum is populated significantly by fermionic models below this scale, and hence the latter serves as our cutoff fc max , 
appearing in (|50| . (fBTjl . A similar situation characterised the ((l+l)-dimensional) boson case in [7J]. The reader is 
invited to compare the order of magnitude of (f6"2"|) with that of (j35l) . 
It is easy to show that 



, 771' 



{8m) 2 



f {0 |: Too(V>i,y>2) :| 0)/ « sin 2 9 — y "''" J AC(7,)I(fc ma x) (63) 
where 

I(fc max ) - /"'"" V - P " (64) 
Jo (p 2 + m^)2 

and we have considered^ <5r7i <C rn. 

Similarly, f (0 \: Tufa, fa) ■ I 0)/ is given by 

f (0 \: Tufa fa) :| 0)/ « - sin 2 fl g-ff^- AC(t ? ) J(fc max ) (65) 

J(fcmax) = / ^ — J-T- (66) 

Jo (p 2 +m^)2 
Representing fc max in units of the characteristic neutrino mass scale m, i. e. 

fcmax = , K > , (67) 



where 



it is easy to show that 



I(kM) = log( K + y/l + k 2 ) 



k(3 + 4k 2 



3(1 +K 2 )2 

t/ — \ , , , /T— -v «(15 + 35k 2 + 23« 4 ) 

J (Km) = log(« + v 1 + k 2 ) ; g 1 . (68) 

15(1 + K 2 )2 



Because of (|56|) . we may write (J63J) as: 



,(0 |:Too(&,&) :| 0) f « sin 2 0^l!^|^C(7 ? )J(fc max ) 



(69) 



The reader is invited to compare this expression with the corresponding one for the boson case, ([56]) , upon taking (p)6"|) 
into account. Upon making the simplifying assumption that the foam is responsible for the whole of the experimentally 
observed mass differences of light neutrino species, we observe that the factor multiplying sin 2 #I(fc max ) is of order 

( (5m 2 ) 3 

•/ C(rj). For late eras, C(rj) ~ 1 (in units of the present epoch scale factor of the Universe). For the biggest 

of the mass differences observed today [20|, Srn 2 ,^ ~ 0.0027 eV 2 (in conventional notation), this factor is of order 
7 x 10~ 118 Mp. Moreover the observed mixing sin 2 #23 contributes factors slightly smaller than one (the current 
experimental data [29j indicate sin 2 (2# 2 3) > 0.87 at 68% confidence level), and, for the ranges of cutoffs considered 
above, the cutoff factors I(fc max ) are of order O(10) at most. The accepted magnitude of the vacuum energy, claimed 
to have been observed today in the form of a positive cosmological constant, A, is A ~ 10~ 122 Mp (Mp ~ 10 19 GeV). 
In order to reproduce such a value one needs 

AtT7 2 

• 2/1 foam -. . ir >— 4 

sin 6 — zih> — ~ 1-4 x 10 . 
m 

This is compatible with other phenomenological tests of space-time foam using neutrinos [lj|. The reader is also 
invited to compare this result with the bound (|28p . derived from stringy uncertainty considerations, in the case of a 
string mass scale M s — O(TeV), upon taking proper account of the theoretical uncertainties due to model dependence, 
as discussed there. 
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The equation of state of this fermionic fluid in the flavour vacuum is approximately determined by: 



WF = - 1 -- ^Hn^ (70) 



(p 2 +mf)2 



From (|70[) . to leading order in dm, we can deduce that wf lies in the range — 1/3<wf<0 since 

±(l ^ , 

3 V 5 3 ( K )(1 + K 2 )i 



where g(n) is a non- negative function given by 

5 ( K )=log( K +Vl + « 2 )--^ r-f (72) 

3(1 + K 2 ) 2 

The asymptotic approach of Wf to — ^ as k — > oo is logarithmic with the cut-off and can be shown to be 

I L 1 fl 7 1 

log(re) t 5 + 8^ ~ 2T 2 



^ = - o ( 1 - TTZTTT ll + O-^H' (73) 



For the cut off of (|62p . ioj? ~ —.17. As k is changed from y2 the value of u>i? rapidly asymptotes to — |, e.g. with 
k = \/20, w F ~ -.27. 

In 3 + 1-dimcnsions, for the fermionic case, the relative factor of 3 in the denominator of the spatial components of 
the stress tensor (as compared to the temporal component (fSTjl ) is due exclusively to the spatial dimensionality. 
In the bosonic case, by contrast, the dimensionality of space is not relevant when evaluating the equation of state 
Wb in terms of the appropriately normal ordered components of the corresponding stress tensor. Hence, in the latter 
case, we obtain % ~ —1, for late eras. 



VIII. DISCUSSION AND OUTLOOK 



In this work we have evaluated above the contributions from (3+l)-dimensional low-energy fermions to the flavour- 
Fock-space vacuum energy and pressure on a brane world punctured by D-particle defects. We have found that the 
pertinent liquid is not describing a cosmological constant vacuum. The equation of state lies in the range > w > — ^ . 
This does not lead to acceleration of the Universe, which requires w < — |. 

Of course, in a D-brane setting, as that of fig. Q] examined here, there are many other contributions to the D3-brane 
world, some of which are notably attributed to bulk D-particles [21j . At late eras (relative to the time of the cosmically 
catastrophic brane collision corresponding to a Big Bang), there are also bulk contributions to the brane dark energy, 
which are due to strings stretched between D-particle defects and the D3 brane world. The reader should recall [21| 
that the contributions to the vacuum energy from the above processes are due to the perpendicular components of the 
relative velocity of the D-particle with respect to the D3-brane world in the bulk; the motion of a D-particle on the D3 
brane world, parallel to its uncompactified components does not lead to any contribution. Hence it is plausible that 
these contributions arc subdominant since in late epochs the motion of the D3 brane in the bulk is extremely slow. It 
is in this sense that the dominant contributions to the D3-brane vacuum energy could come from the above-described 
processes of capturing and splitting of flavoured string states on the brane, corresponding to the flavour Fock-space 
vacuum contributions evaluated above. 

To ensure that the flavour vacuum contributions to the dark energy leads to accelerating Universes at late epochs, 
as the current phenomenology indicates, one should have bosons simultaneously present with fermions. In fact the 
bosonic and fermionic contributions to the vacuum energy and pressure are algebraically additive. In the case of 
D-particle foam, only electrically neutral particles interact non trivially with the D-particle defects. 

In the context of supersymmetric low-energy field theories, such as those derived in the low-energy limit of super- 
strings, the relevant bosons may be the sneutrinos, the supersymmetric partners of neutrinos, which have large masses 
due to target-space supersymmetry breaking. However, the relative mass differences between mass eigenstates may 
be assumed sufficiently small, since the mass differences are independent of supersymmetry, especially if, according to 
our D-particle foam model, they are quantum-gravitational in origin. Hence even if the partners have a much greater 
mass, due to supersymmetry breaking, we may assume that, among different flavours, the same small mass differences 
that characterize the fermionic excitations also characterize the bosonic superpartner flavours. In this sense, one has 
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contributions to the vacuum energy density and pressure from the bosons, which are of the same order as those of 
fermions. In realistic supersymmetric models, the total equation of state may be complicated, as it depends on the 
various fluids that participate in the flavour vacuum structure. Nevertheless, it is possible to have an equation of 
state that guarantees a late-era acceleration of the universe. 

However, technically the extrapolation of the above results to supersymmetric cases is not a trivial task. A super- 
symmetric theory is by construction typically an interacting theory, while above the quantization procedure adopted 
was based on free excitations. We hope to come back to this important issue in a forthcoming publication. 
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